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ABSTRACT 

We propose that giant flares on Soft Gamma-Ray Repeaters produce rela- 
tivistic, strongly magnetized, weakly baryon loaded magnetic clouds, somewhat 
analogous to solar coronal mass ejection (CME) events. Flares are driven by 
unwinding of internal non-potential magnetic fields which leads to slow build-up 
of magnetic energy outside of the neutron star. For large magnetospheric cur- 
rents, corresponding to a large twist of external magnetic field, magnetosphere 
becomes dynamically unstable on Alfven crossing times scale of inner magne- 
tosphere, tA ~ Rns/c ~ 30/isec. The dynamic instability leads to formation 
of dissipative current sheets through development of tearing mode (Lyutikov 
2003). Released magnetic energy results in formation of a strongly magnetized, 
pair-loaded, quasi-spherically expanding flux rope, topologically connected by 
magnetic field to the neutron star during the prompt flare emission. Expansion 
reaches large Lorentz factors, F ~ 10 — 20 at distances r ~ 1 — 2 x 10'' cm, where 
lepto-photonic load is lost. Beyond this radius plasma is strongly dominated by 
magnetic field, though some baryon loading, with M <^ E/c^, by ablated neu- 
tron star material may occur. Magnetic stresses of the tied flux rope lead to 
late collimation of the expansion, on time scales longer than giant flare duration. 
Relativistic bulk motion of the expanding magnetic cloud, directed at an angle 
6 ~ 135° to the line of sight (away from the observer), results in a strongly non- 
spherical forward shock with observed non-relativistic apparent expansion and 
bulk motion velocities jSapp ~ cot^^/2 ~ 0.4 at times of first radio observations 
approximately one week after the burst. Interaction with a shell of wind-shocked 
ISM and then with the unshocked ISM leads to deceleration to non-relativistic 
velocities approximately one month after the flare. 
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1. Introduction 

Magnetar emission {see, e.g. , Woods & Thompson 2004, for review) is powered by 
dissipation of a non-potential (current-carrying) magnetic field (Thompson et al. 2002). 
Dynamo mechanism operating at birth of neutron stars creates a tangled magnetic field 
inside a neutron star, which is prevented from unwinding by rigidity of the crust. Current- 
carrying plasma exerts Lorentz force on the crust, which is mostly balanced by lattice strain. 
For strong enough magnetic fields, Lorentz force may induce a stress that exceeds the critical 
stress of the lattice. This leads to crustal motion, which should occur along cquipotcntial 
surfaces: crust will be rotating. Crustal rotation and associated twist of magnetic field 
lead to expulsion of the electric current from inside of neutron star into magnetosphere. 
Dissipation of magnetospheric currents is responsible for persistent emission (Thompson et 
al. 2002), while sudden reconfiguration of magnetic field may produce fiares (Lyutikov 2003). 
Giant flare of SCR 1806-20 on December 27 2004 (we will refer to it as "the GF") puts new 
constraints on the model that we discuss in this paper. 

2. Where was energy stored right before the flare? — In the magnetosphere 

2.1. Short rise time 

GFs are powered by dissipation of magnetic field energy. One of the principal issues is 
where most of the magnetic energy has been stored prior to the GF: in the magnetosphere 
or in the neutron star crust. These two possibilities are related to two models of giant 
flares of the SGRs. First, a giant flare may result from a sudden untwisting of the internal 
magnetic field Thompson & Duncan (1995, 2001), TD95 and TDOl below. In this case, 
a large and quick (on time scale on a fiare) rotational displacement of the crust leads to 
increased twisting of magnetospheric magnetic field lines. Alternatively, slow untwisting of 
the internal magnetic field leads to gradual twisting of magnetospheric field lines, on time 
scales much longer than GF, until it reaches a dynamical stability threshold due to increasing 
energy associated with current-carrying magnetic field. Then sudden relaxation of the twist 
outside the star and associated dissipation and magnetic topology change lead to flares, in 
analogy with Solar flares and Coronal Mass Ejections (CMEs). Note, that even in case of 
crustal storage of magnetic field energy before the flare (TD95, 2001), dissipation also occurs 
in the magnetosphere, not in the crust. 

The best test of these two alternatives is the time scale for the development of a flare. 
Since energy involved in the GF requires that a large fraction of the magnetosphere is 
affected, a typical size of the active region is of the order of the neutron star radius. Sudden 
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unwinding of the crust should occur either on shear wave or Alfven wave crossing time of the 
star tA,NS ~ Rns/Vs,ns ~ Rns/yA^NS ^ 0.2—0.5 s (TD95), while magnetospheric instability 
may develop on time scales as short as Alfven crossing time of the inner magnetosphere, 
tA,ms ~ Rns/c ~ 30/US (Lyutikov 2003). Observations of the December 27 GF show very 
short rise time ~ 0.25 msec (Palmer et al. 2005). Very short rise time of the GF points to 
magnetospheric origin of GF (in a sense that right before the flare the energy to be released 
is stored in the magnetosphere). 

2.2. Pre- and Post-burst evolution of persistent emission 

XMM-Newton observation of SGR 1806-20 months before the GF shows an increased 
activity, with persistent flux increasing by a factor of two, spectrum hardening (photon power 
law index decreased from 2.2 to 1.5) and spindown increasing (Mereghetti et al. 2005), all 
in agreement with the prediction of twisted magnetosphere model (Thompson et al. 2002), 
implying an increasing twist before the GF. In addition, two months after the GF pulsed 
fraction and the spin-down rate have significantly decreased and the spectrum softened (Rea 
et al. 2005). The same occurred in SGR 1900+14 following the August 27 GF (Woods et al. 
2001). All of these effects are consistent with increasing of the twist during the time leading 
to the GF and decreased twist of external magnetic fields after the GF, brought about by 
reconnection: (i) in the rcconncction model the post-fiare magnetosphere is expected to have 
a simpler structure, as the pre-fiare network of currents has been largely dissipated; (ii) non- 
thermality of the spectrum is a measure of the current strength in the bulk of magnetosphere 
(Thompson et al. 2002; Lyutikov & Gavriil 2005) with softer spectra corresponding to a 
smaller twist; (iii) spin-down rate depends on the amount of current fiowing through open 
field lines and is smaller for a smaller twist. Note, that since open field lines occupy only a 
small fraction of magnetosphere, spin-down rate probes current in a relatively small region 
which should, on one hand, correlate with typical current in the magnetosphere on long time 
scales, but on the other hand may show large deviations on short time scales. 

2.3. Ejecta must carry a lot of magnetic field 

Magnetospheric storage and release of energy leads to the following consequences. First, 
it is hard to see how most of the energy released in the magnetosphere can be spent on 
heating the surface of the neutron star and generating heavy ion-loaded outflows. Secondly, 
dissipation of magnetic field cannot create magnetic field-free plasma: it is likely to be 
limited to equipartition fields since at this point induced magnetic field of gyrating relativistic 
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particles will create magnetic field comparable to the initial field (at temperatures near 500 
keV energy in photons will be comparable to energy in electrons) . Thus, only approximately 
half of the magnetic field is expected to be converted to particle and photon energy, not 
much more. Most of the dissipated energy will be later radiated away and/or spent on pdV 
work during expansion. 

Thus, magnctospheric release of energy, indicated by very short rise time of the GF, 
decreasing persistent emission and softer post-fiare spectrum, leads to the conclusion that 
expanding plasma must be strongly magnetically dominated. In this paper we examine 
consequences and consistency of the model based on these premises. 

3. Overview of the model 

Before discussing various details, let us present a short overview of the model, which 
qualitatively resembles models of Solar flares with the difference that Sun supports an actively 
operating dynamo, while in magnetars dynamo operated only during birth of a neutron 
star. The energy that will be released in a GF is initially (at times long before the flare) 
stored in electric currents flowing inside neutron star. These currents arc generated during 
birth of the neutron star and are slowly pushed out into magnetosphere, gated by slow, 
plastic deformations of neutron star crust. This creates active magnctospheric regions, in 
analogy with Solar spots. An active region consists of a sheared arcade of magnetic flux 
and surrounding non-potential magnetic structures. As currents (and with them magnetic 
energy and helicity) are pushed outside the neutron star, magnetosphere adjusts slowly to 
the changing boundary conditions. During this phase magnetic energy is slowly stored in 
the magnetosphere. As more current is pushed outside, magnetosphere reaches a point of 
dynamical instability beyond which stable equilibrium cannot be maintained. Grossing the 
instabilities threshold leads to changing of magnetic configuration on Alfven crossing time 
scale, formation of narrow current sheets, and onset of magnetic dissipation (this process 
is sometimes called magnetic detonation {e.g. Gowley and Artun 1997)). This has two 
consequences. First, a large amount of magnetic energy is converted into kinetic plasma 
energy and photons. Secondly, dissipation allows to change magnetic topology and leads to 
formation of an expanding magnetic loop that eventually break away from the star. Initially, 
after onset of reconnection, kinetic pressure of optically thick pair plasma and magnetic 
stresses are comparable, so that expansion is quasi-isotropic and reaches rclativistic Lorentz 
factors > 10, determined either by baryon loading or, in case of very small baryon loading, 
by the amount of residual pairs. 

During the prompt phase of the GF, the expanding magnetic loop remains attached to 
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the star, see Fig. 1. Tying of the expanding loop to the star eventually leads to coUimation 
of the explosion into a wide opening angle, of the order of one steradian. After losing the pair 
load, the expanding cloud is dominated by magnetic field (magnetic cloud) and eventually 
disconnects form the neutron star, moving relativistically away from the observer at an angle 
~ 135° degrees. This results in apparent subluminal proper expansion and proper velocity. 
Eventually, energy of the magnetic cloud is transfered to the strongly anisotropic forward 
shock which produces the observed afterglow radio emission. 



4. Prompt and tail emission of the GF 

4.1. Twisting of external magnetic field and GF precursors 

The central point of our suggestion is that winding-up of external magnetic field pro- 
ceeds on long time scale, much longer than GF. The winding can occur on time scale of 
approximately two months before the GF, indicated by increased activity of the source 
(Mereghetti et al. 2005). Alternatively, onset of fairly rapid plastic deformation of the crust 
may be related to weak emission events (precursors) that seem to precede GFs. [In case of 
GF of SGR 1900+14 a precursor was seen ~ 0.4 seconds before the main burst (Feroci et al. 
2001) while in case of SGR 1806 a relatively powerful event occurs approximately hundred 
fourty seconds before the main burst (Palmer et al. 2005).] During the quiescent period 
between the precursor and GF, a patch of a crust is continuously rotated by the Lorentz 
force, balanced both by elastic and viscous stresses in the crust. These two possibihties 
are not mutually exclusive: slow evolution on time scale of months may be followed by a 
relatively fast twisting over hundreds of seconds before the GF. 

Consider a crustal plate of size R rotating under the influence of Lorentz force, balanced 
by viscous stresses at the base of the crust. The dissipated power is (Landau & Lifshitz 1975) 

... ^ 1^1^ ^ .3 . (^) ' (^) [^f ,1) 

/ N 5/4 

where v ~ 10"^ ( lo^^^xm--^ ) viscosity of neutron star (Cutler & Lindblom 1987), p 

is density at the base of the crust, T^ot = 140 (^) s is the rotation period of the plate, 
taking into account that instability occurs after rotation of A0 radians. Total viscously 
dissipated energy is Eyi^ ~ L^iscT ~ 2 x 10^^ erg. [For SGR 1900-^14 with Trot = 0.4 sec 
Lvisc ~ 3 X lO^^ergis"^ and E'^s ~ 10^^ erg, but a very short T^ot may indicate that twist 
was near critical before onset of rotation, A0 <^ 1.] Since this energy is released deep in 
the crust, where thermal diffusion time to the surface is much longer than T^ot, most of the 
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heat is absorbed by the core (c/. Lyubarsky et al. 2002) and does not show as increased 
persistent emission between the precursor and the main flare. 

During plastic creep, elastic strain is much larger than plastic strain. The difference 
between the two can be expressed in terms of how much magnetic fleld exceeds fleld at the 
critical strain: 

AB ^ 2^>V>^VH,.^"T;r - 2.7 X 10'»G {j^f {^f (2) 

{AB = 2.2 X 10^2 G for SGR 1900+14). Thus, during plastic creep only a small fraction of 
magnetic energy is dissipated in the crust, most of it is pumped outside of the star. 

Note that in case of the Sun, TRACE satellite has detected rotation of sunspots asso- 
ciated with largest (X-class) flares days before the flare (Mewaldt et al. 2005). In addition. 
Solar CMEs also starts before the accompanying X-ray flare: there is a quiet growth period 
of approximately 30 minutes before formation of the dissipative current sheet (Zirin 1988). 



4.2. During the giant flctre expansion must be relativistic 

Relativistic expansion at the time of the GF follows from the conventional compactness 
argument. For luminosity L ~ 10^^ erg/s and variability time scale ~ 1msec, optical depth 
to pair-production is (see also Nakar et al. 2005) 

-^7-7 ~ . ^^^r. 2 X 10^^ (3) 

If plasma were to remain non-relativistic, photon diffusion times would be long, inconsistent 
with short observed variability time scale. This estimate immediately excludes large baryon 
loading: M must be ^ E/(?. This has always been a standard view of giant flares, (Feroci 
et al. 2001; Thompson & Duncan 2001). [A possibility that initial 7-ray spike was produced 
in a relativistic outflow without dynamically important magnetic field and did not contribute 
significant amount of energy to afterglow is hard to reconcile with magnetospheric origin.] 

At early stages dissipation of magnetic field creates an optically dense lepto-photonic 
plasma emerged in magnetic field with T ~ (L/ ATrRf^gasrY^'^ ~ 300 keV (see also Nakar 
et al. 2005). Qualitatively, quasi-spherical expansion of a strongly magnetized pair bub- 
ble resemble the unmagnetized case, but there are important differences in the asymptotic 
dynamics, outlined in Appendix A. Initially, plasma expands with bulk Lorentz factor in- 
creasing approximately linearly with radius F oc r, while rest temperature decreases T ~ 1/r. 
After reaching T-t ~ 20 keV at r± ~ 1.5 x 10^ cm (at which point F ~ 15), plasma becomes 
optically thin. Observed emission is thermal with To^g — TT^ ~ 300 keV. 
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The main implication of relativistic expansion is that outflow is not heavily loaded with 
baryons. In §6 we will show that this picture is consistent with observed non-relativistic 
expansion velocities of the afterglow. 

4.3. Initial millisecond spike is nearly isotropic 

If initial spike were strongly anisotropic, with luminosities inside and outside some 
emission cone different by orders of magnitude, then we would have observed many more 
tails without initial spikes. Not a single tail without a spike has been seen. Tails are obviously 
only weakly anisotropic, emitted by a trapped fireball (TD95), while their intensity is well 
above threshold of detectors. 

Since initial 7-ray emission is nearly isotropic, it is unlikely to be produced by a strongly 
jetted outflow (contrary to Yamazaki et al. 2005). It is still feasible that the initial spike is 
weakly anisotropic, with radiation intensity and Lorentz factors changing by some factor < 2 
depending on direction. [Note, that tail emission in all cases was of the same order, while 
the energies of the initial spikes were vastly different. This fact is, on one hand, consistent 
with some structured jetted emission of the initial spike, so that all bursts are the same but 
initial spike is viewed from different angles (Yamazaki et al. 2005), but on the other hand 
it contradicts the fact that afterglow emission in case of the GF was several orders brighter, 
arguing in favor of larger total energetics. Constant tail emission may be explained as a 
limiting effect of magnetar magnetic field: above some threshold flare energy, the amount of 
trapped plasma depends only on the strength of confining poloidal magnetic field and not 
on the amount of the released twist.] 

4.4. Time scales of the GF 

There are several time scales in the model: first is the slow initial twist of external 
magnetic field §4.1. Secondly, there are several time scales associated with the GF itself: 
(i) sub-millisecond initial rise .25 msec (ii) ~ 5 msec rise to the main peak, (iii) hundreds 
millisecond total duration of the spike, (iv) tens of seconds tail emission. Finally, there 
is afterglow time-scale from one week to ~ one month, the latter we identify with non- 
relativistic transition of the ISM blast wave. 

Let us discuss the time scales of the GF itself. Primarily, we associate the shortest 
time scale observed in the burst ~ 0.25 msec with the Alfven crossing time of the in- 
ner magnetosphere. It reflects the dynamical evolution of the magnetosphere after it has 
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crossed stability threshold. Thus, the very first photons are emitted while plasma is still not 
expanding relativistically. 

Secondly, since emission requires dissipation of energy, dissipative time scales become 
important as well. In magnetar magnetospheres development of dissipative tearing instability 
occurs on time scale intermediate between Alfven times scale and resistive time scale: 
ttearing ~ \/vc/Rns whcrc T] IS plasma resistivity. If resistivity is related to plasma skin 
depth, f] ~ c^/ujp, where Up is plasma frequency, the growth rate of tearing mode is ~ 10 
msec for a current sheet of width ~ -Rats (Lyutikov 2003). Intermediate time scale ~ 5 msec 
observed in the GF may be related to development of tearing mode in the current sheets 
formed during onset of dynamical instability. 

We associate the overall duration of the spike ~ 100 msec with the dynamical time of 
the expanding magnetic cloud ~ 2T'^c/Rns ~ 25msec for F ~ 20. This is the minimum 
time it takes for expanding strongly magnetized bubble to come into causal contact with 
itself. One expects that on this time scale magnetic cloud re-adjusts its internal structure 
and relaxes to a minimal energy state. This relaxation occurs, e.g. , through reconnection, 
which in relativistic case may proceed with the inflow velocity reaching the velocity of light 
(Lyutikov & Uzdensky 2003). 

Finally, typical time scale for the evolution of the tail of the GF, tens of seconds, is well 
described by radiation leaking from a plasma trapped on closed magnetic field hues, TD95. 

4.5. Quasi-thermal spectrum of the initial spike. 

Thermal spectrum results from radiation escaping from a (strongly magnetized) fireball 
that becomes optically transparent. Observer temperature of hundreds of keV corresponds to 
rest frame temperature of ~ 20 keV, when plasma becomes optically thin to pair production, 
boosted by Lorentz factor ~ 10 (Goodman (1986); Paczynski (1986), see also Nakar et al. 
(2005)). 

4.6. Mass loading and terminal Lorentz factor 

In appendix A we consider dynamics of a hot, strongly magnetized expanding flow 
carrying toroidal magnetic field. Flow is accelerated by magnetic and pressure forces, while 
both matter inertia and magnetic field energy density provide effective loading of the fiow. 
In addition, in case of a large scale magnetic field considered here, there is extra conserved 
quantity: magnetic fiux. This plays an important role in the overall dynamics of the fiow 
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(c/. Kennel & Coroniti 1984). 

If source luminosity is L, mass loss rate is Mq and EMF is £ (these are conserved 
quantities), then the terminal Lorentz factor and terminal magnetization parameter are 

- ^ - 

Which, formally, expresses the fact that magnetic field provides additional effective loading 
(factor 1 + (Too), but the amount of loading depends non-trivially on parameters of the 
flow. In the strongly relativistic limit, outflow typically reaches Alfven velocity (in fact 
fast magnetosonic), at which point Lorentz factor is related to the terminal magnetization 
parameter as (Michel 1971; Goldreich & Juhan 1969) 

Too = V Coo) (5) 

To estimate the maximum possible Lorentz factor, we note that the minimum mass 
loading is determined by residual left-over pairs, determined by equating annihilation and 
expansion rates (Goodman 1986; Paczynski 1986; Nakar et al. 2005) 



If most energy is in magnetic form, this corresponds to (Tmax = E/M„iinC^ = 4.5 x 10^. The 
corresponding maximum Lorentz factor is Vmax = \/<^max = 6.7 x 10^. Lower limit on F 
comes from the observed thermal temperature of the initial spike: 

Tmin ~ ^ ~ 10 - 20 (7) 

-tit 

Thus, the flow must be only weakly polluted by baryons. This picture is in full agreement 
with TDOl. In what follows we adopt a minimum value of F = 10 for numerical estimates. 
Then, estimating E ~ Lts and M ~ Mtg and using Eq. (5), Eq. (5) gives 

M-j§^,-^-l0^h (8) 
This is the upper limit on amount of mass ejected during the GF. 



4.7. Plasma physics issues 



The proposed model of the GF builds on the models of Solar CME and, similarly, 
has a number of problematic plasma physics issues, (see, e.g. Priest & Forbes 2002, for 
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review) One is what is known as Aly-Sturrock paradox (Aly 1984; Sturrock 1991): opening 
of field lines, which is necessary to generate an outflow, requires an increase in the magnetic 
energy in the system, while the storage model of CMEs requires magnetic energy to decrease 
during formation of magnetic cloud. There is a number of ways Aly-Sturrock paradox can 
be circumvented, most important being magnetic reconnection which can change topology 
of the field line. 

Another problem is that injection of current occurs on a finite amount of magnetic flux 
so that in order to expand the newly formed magnetic cloud has to break through overlying 
closed dipolar flcld lines. This is achieved by reconnection at the null point at the leading 
edge of the magnetic cloud, cf. "magnetic breakout" model of Antiochos et al. (1999). 
Reconnection transfers unshared magnetic flux associated with overlying dipolar fleld to 
neighboring flux tubes, allowing the sheared fllament to expand and erupt outward. The 
rate of reconnecting adjusts so that radial (as seen from the star) propagation velocity is the 
Alfven velocity, which in this case is nearly the velocity of hght. 

5. Radio afterglow: qualitative description 

5.1. Expansion in SGR wind 

As magnetic cloud expands, it becomes transparent at r± and its pair density falls by 
many orders of magnitude. At this point magnetic cloud becomes strongly magnetically 
dominated. Initially magnetic cloud is topologically connected to the star, but eventually 
reconnection should happen at the footpoints of the magnetic field lines, disconnecting the 
magnetic cloud from the star. At this point magnetic cloud d starts to expand into preexisting 
SGR wind. It is expected that SGR wind is strongly relativistic, with Lorentz factors 
» 10 — 20. Thus, for intermediate baryon loading (such that magnetic cloud expands with 
r ~ 10 — 20) magnetic cloud never overtakes the wind, so that expansion occurs as if in 
vacuum. Most of the magnetic energy is concentrated in a shell close to the bubble surface 
with thickness of the order of cts ~ 10^ cm, where tg ~ 100 msec is flare duration. 

5.2. Apparent constant non-relativistic expansion velocity is due to 
relativistic strongly anisotropic expansion. 

Observations of constant expansion velocity from two to flve weeks after the burst have 
been interpreted as evidence in favor of large baryon loading, and, as a consequence, weak 
relativistic initial expansion velocities (Granot et al. 2005). A we argued above, this cannot 
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be the case due to compactness constraints: flow must be relativistic with small baryon 
loading. 

Apparent non-relativistic expansion velocity can be due to relativistic anisotropic ex- 
pansion with little emission within the cone l/F to the line of sight. If emitting mate- 
rial is moving relativistically at an angle 9 ^ l/F, then apparent expanding velocity is 

Papp = P sin 9/{l — (3 cos ^) ~ /? cot 9/2. To reproduce observed Papp ~ 0.3 — 0.4, it required 
that outflow is directed away from the observer at an angle ~ 135°. 

We have arrived at a seeming contradictory picture: initially, during the 0.2 sec spike 
of the GF, expansion should be nearly isotropic, while at later times, > 1 week expansion 
is strongly anisotropic. Thus, magnetic cloud should become strongly anisotropic between 
0.2 seconds and 7 days. It is unlikely that anisotropy is achieved by internal magnetic 
stresses of the freely expanding magnetic cloud: in relativistic case self-collimation is strongly 
suppressed {e.g. Bogovalov 2001). Similarly, as we argue in §5.3, coUimating effects of the 
dipolar magnetic field cannot account for strong anisotropy given the nearly spherically 
symmetric initial explosion. We propose that expanding magnetic cloud becomes strongly 
anisotropic due to fact that magnetic fields of the cloud remain attached to the neutron 
star during most of the prompt phase. Thus, at this times the magnetic topology of the 
expanding plasma is that of a fiux rope, see Fig. 1. 



5.3. Spheromac or flux rope? 

Last decade two models were proposed for the structure of interplanetary magnetic 
clouds ejected from the Sun: a magnetic fiux rope and spheromac. The principal difference 
between the two is that in case of spheromac, magnetic cloud disconnects from the magnetic 
field of the Sun at early stages of ejection and becomes quasi-spherical {e.g. Gibson & Low 
1998), while flux rope remains attached to the Sun for a very long time (even at the orbit 
of the Earth). The two models lead to very different dynamics of the magnetic clouds. 
Presently, spheromac model seems to be inconsistent with data {e.g. Farrugia et al. 1995), 
while magnetic ffux rope model explains well the internal magnetic structure ejected into 
interplanetary space (Marubashi 2000). 

In case of the GF, spheromac model seems to be inconsistent with data for the following 
reason. As we argued above, the initial explosion should be quasi-isotropic, while at later 
times it should become strongly anisotropic and one-sided. It unlikely that relativistic, 
strongly anisotropic explosions are produced due to coUimating effects of internal magnetic 
field of the expanding blob which is disconnected from the star: for relativistic expansion 
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coUimation by internal magnetic field is kinematically suppressed {e.g. Bogovalov 2001), 
and in any case cannot produce one-sided explosion. Spheromac also cannot be efficiently 
collimated by external dipolar magnetic field, since inside spheromac internal kinetic and 
magnetic pressures scale as oc -Bfp/^ ~ {Bgph is a typical magnetic field inside spheromac), 
while -Bjjpo/ar ~ Thus, if the plasma to be ejected disconnects from the stellar magnetic 
field early on, one may expect only weak coUimation and, as a result, weak relativistic bulk 
motion. Since the overall expansion should be strongly relativistic, as we argued in §4.2, 
observed emission will be dominated by parts of the shock moving towards the observer, 
with somewhat smaller Lorentz factor than average, but having large apparent velocity. To 
illustrate the point, in Fig 2 we plot Doppler factor 5 — l/(r(l — Pcos9) and apparent 
transverse velocity for a relativistic shock expanding with Lorentz factor F = 3 and moving 
at Sob = 135° with bulk Lorentz factor F;,^;^ = 2 as a function of the angle 6 between the 
explosion direction and emission point in the rest frame of the explosion and emission points 
located in the plane containing direction of bulk motion and line of sight. Points on the shell 
with highest Doppler boosting have large apparent transverse velocities, in contradiction 
with observations. 

On the other hand, expanding magnetic cloud confined by a fiux rope (with B^^,^^ ~~ 
r""^) may in principal provide coUimation (this would correspond to only regions near 9ob < tt 
in Fig. 2 contributing to observed emission). Details of this late coUimation need to be 
investigated numerically. From observational point of view, at late times expansion should 
be confined to a fairly broad angle, of the order of tt radians, and not to thin, GRB-like jet. 



6. Bcisic afterglow parameters 
6.1. Geometry 

For numerical estimates we chose Ooj, — Sn/A — 135°. Then apparent velocity Papp ~ 
BcotO/2 = 0.41/5, Doppler factor S ~ l/(2rsin^ e/2) = 0.58/F and observer time IS given 
by T = 2rsvD?9/2/{f3c) = 1.70r/(/9c). We assume a strongly magnetized fiow with total 
isotropic energy Ef.j ~ 10"^^ erg which reaches terminal Lorentz factor F = 10, and is colli- 
mated into angle (ii7/(47r) ~ .1 (so that "typical opening angle is ~ 36°). We also normalize 
surrounding density to n = 1 cm^. 
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6.2. Typical radii and time scales 

Before the explosion, the magnetar is surrounded by a nearly empty bubble blown by 
the magnetar wind with a typical size 

r. ^ f--^V'' ^ 1.2 X 10-cm f^)"'' V' (9) 

\A-KnmpcV^J Vlcm-s/ VlOOkms^V 

where ~ 10^"^ erg/s is average spindown luminosity of SGR 1806 and Vns is velocity of 
the neutron star. Since magnetar wind is expected to have Lorentz factor ^ 10, until 
magnetic cloud expands freely, without slowing down. Note, that Vg is the minimum distance 
from the neutron star to the shell, depending on the relative orientation of the neutron star 
velocity and direction of explosion the time when magnetic cloud overcomes the shell can be 
larger by factor of ~ 2. 

Were it to expand in a constant density medium, the flow would starts to decelerate at 

Since r^ec < f^s, magnetic cloud starts interacting with the a shell of shocked ISM plasma at 
r ~ r,. In observer time this occurs at 



s- 



T ~ 1.7-^ = 8.2days (11) 

Since the amount of the mass contained in a shell of shocked ISM plasma is ~ nniprg, 
after encountering the shell Lorentz factor of the magnetic cloud falls down to 

1/2 / P \ 1/2 / T/-..„ \ 3/2 , ^ , 1/4 / ^ 



( y'_^^j E Y ( V ( ^ V 

\dQnmpC^r^ J ' yw^^ergs'^ J VlOOkms~V Vlcm-^/ 



0.1 X Att 



(12) 



Doppler beaming factor at this point is 5 = 0.2. 

Transition to the non-relativistic expansion occurs at 

/ 3£' \ 1/3 

■■"'-(^S^) -l-^xlO"- (13) 
corresponding to the observer time 

Tnr = 16.5days (14) 

Tnr is an estimate of time when the velocity of the blast wave starts to deviate considerably 
from c. It typically takes two times longer for the velocity to fall below 0.5 c, and for the 
Doppler factor to become within 15% of unity. One expects to see a peak of emission at 
the moment when Doppler de-boosting becomes insignificant, approximately at ~ 2T„r ~ 33 
days. 
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6.3. Total energy 

Form the standard cquipartition argument (Pacholczyk 1969), and taking into account 
Lorentz transformation of flux, oc and frequency, oc 6, the minimal energy of relativistic 

electrons plus magnetic field is 

where a ~ 0.5 is the spectral index. This estimate is by a factor 5-2+40/7 ^ larger than 
the one based on assumption of non-relativistic expansion. Corresponding magnetic field (in 
laboratory frame) is 

..„^o.c(^)--.-(^)^"(^)-(i;^)"'"^"". (-) 

which is larger by a factor (5-i-2a/7 ^ g than the one based on assumption of non-relativistic 
expansion. Note, that minimal energy argument address only energy in magnetic field and 
relativistic electrons. It is expected that most energy of the forward shock resides in protons, 
so that the total energy in the outflow may be order(s) of magnitude larger than (15), bringing 
it in line with the total energy released in 7-rays. 



7. Discussion 

In this paper we outlined a model of magnetar GF based on the analogy with Solar 
Coronal Mass Ejections. Very short rise time scales of the GF indicate that GF flare is 
driven by dissipation of energy stored in the magnctosphere right before the burst (as sug- 
gested by (Lyutikov 2003)) and not in the crust of the neutron star (as proposed by TD95). 
Initially, the explosion is loaded with pairs and is quasi-isotropic. Magnetic field topology of 
expanding plasma resemble flux rope model of CMEs, which leads to late time coUimation 
and anisotropic expansion. Expanding magnetic cloud is strongly dominated by magnetic 
field, weakly baryon loaded, with M <^ E/c^, and strongly relativistic. Weeks after the flare, 
magnetic cloud still expands relativistically, strongly anisotropically and is moving away from 
the observer, resulting in apparent expansion velocity Papp ~ 0.4. At approximately 33 days, 
corresponding to a bump in the light curve, expansion velocity falls below c/2. The main 
prediction of the model is that we may see a medium-energy flare with a very bright radio 
afterglow, when the explosion will be beamed towards the Earth. Statistically, it should take 
approximately ~l/dQ~10GFs. 
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Magnetospheric dissipation following plastic deformation of the crust is also consistent 
with suggestion of Jones (2003) that that neutron-star matter cannot exhibit brittle fracture 
and should instead experience only plastic deformations. On the other hand, complexity of 
earthquakes, especially so called deep focus earthquakes occurring at high pressures, indicates 
that Jones (2003) argument is not the end of the story. For example, crust response may 
depend on value of a strain, being plastic at low strains and brittle at high strains (Frohlich 
1989). 

Can crustal fractures model be consistent with above arguments? One possible way is 
to invoke small scale, ~ 100 meters, initial crustal deformation (so that the rise time of the 
GF is short enough), which triggers larger scale deformations in an avalanche- type process 
(TD95, TDOl). In addition, relatively bright and long-hved afterglow following August 27 
flare is well fitted by the deep crustal heating model (Lyubarsky et al. 2002). The crustal 
fracture model also has a better chance of explaining post-flare activity ( August 27 and 
March 5th flares followed by a burst-active period) in analogy with earthquake aftershocks. 

One possible way to distinguish the models is that reconnection-type events may be ac- 
companied by coherent radio emission resembling solar type-Ill radio bursts Lyutikov (2002). 
The radio emission should have correlated pulse profiles with X-rays, narrow-band-type ra- 
dio spectrum with Au < v with the typical frequency v \ GHz, and a drifting central 
frequency. This requires catching a burst in simultaneous radio and X-ray observations. 

We would hke to thank Roger Blandford, Yuri Lyubarky and Christofer Thompson for 
numerous discussions. 
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A. Dynamics of magnetized pair-loaded flows 

In this section we consider dynamics of a hot relativistic, spherically symmetric, station- 
ary outflow carrying toroidal magnetic field. Since the energy release during GF lasts ~ 0.1 
s, on small time scales and for radii < 10^ cm flow may be considered stationary. In addition, 
internal energy density in the expanding magnetic cloud scales slower with the radius, oc r~^, 
than energy density of the dipolar field oc r~^, so that dynamical effects of dipolar magnetic 
field quickly become negligible. At this stage expansion is quasi-isotropic, as discussed in 
§4.3. Rotation and dynamical effects of the poloidal magnetic field are neglected. 

The dynamics of such a warm magnetized wind is controlled by three parameters: energy 
L, mass flux M and the electro-motive force £ produced by the expanding magnetic cloud. 
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Total energy flux may be divided in two forms: mechanical Lm and electromagnetic Lem 
luminosities. 

L^Lm + Lem (Al) 

We wish to understand how the parameters of a fully relativistic flow (velocity /3, pressure p, 
magnetization cr) evolve for an arbitrary ratio of both Lem/ Lm and p/ p {p is the rest- frame 
mass density). 

The asymptotic evolution of the flow is determined by conserved quantities which may 
be chosen as the total luminosity L, the mass flux M and the EMF £. Thus, the central 
source works both as thruster and as a dynamo. 

The formal treatment of the problem starts with the set of relativistic MHD equations 
which can be written in terms of conservation laws. In coordinate form and assuming a 
stationary, radial, spherically symmetric outflow with toroidal magnetic field, relativistic 
MHD equations give 

^dr [r\w + b^) ^r^]^0 (A2) 

^dr [r' {{w + b')P'r' + {p + 6V2))] - y = (A3) 

^dr [rbpT] = (A4) 

^dr [r'ppr] = (A5) 
The above relations can be simplified by defining 



(A6) 



where we assume that fluid is polytropic with adiabatic index Ta- w = p + y^-zjP, £ is the 
electromotive force. 

It is convenient to introduce two other parameters: the magnetization parameter a as 
the ratio of the rest-frame magnetic and particle energy-density and a fast magnetosonic 
wave phase velocity f3f 

^^^^L^^' ^/ = iT^ + (rT^ = ^^«-^)l^^-^J+^2-^«^L^ 

Using the three conserved quantities L, M and E the evolution equation becomes 

1 (r„ - 1) (l3L - PVM - eA 

-AT = ^ ^ ^ (A8) 



' ^ (^pL{p2 + 1 - r„) + (r„ - i)/3rM - (2 - r^)^:^) 
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Eliminating £ in favor of F/ we get a particularly transparent form for the evolution of 
Lorentz factor 

(£l^a r = 

D^P ' (u.-r.p)r ^ ' 

Equation (A9) is nozzle- type flow (e.g. (Landau & Lifshitz 1975)). The Ihs of eq. (A9) 
contains a familiar critical point at the sonic transition T = Tf. The positively defined 
rhs describes the evolution of Lorentz factors due to kinetic pressure effects. In the case 
of purely radial expansion the magnetic gradient forces are exactly balanced by the hoop 
stresses, so that magnetic field does not contribute to acceleration. Prom Eq. (A9) it follows 
that super-fast-magnetosonic fiows accelerate while sub-fast-magnetosonic fiows decelerate. 
It also follows that terminal Lorentz factor of the flow is determined by the condition drV — 
which implies that either p = or f3oo = 0. Condition j3oo = can be reached only 
for subsonic flows with £ = Q. Neglecting the /9oo = solution, the terminal velocity of 
magnetized flow is determined by the condition: 

L = r^M+|i (AlO) 

Poo 

For each set of parameters (L, M and £) there are generally two solutions of Eq. (AlO) 
for the terminal four-velocity. The only exception is the case of zero magnetization, £ = 0, 
when the terminal (supersonic) Lorentz factor is uniquely determined by Foo = L/M. 

For non-zero magnetization the terminal velocity cannot be determined uniquely from 
given L, M and £. Solutions exist only for L/Af larger than some critical value L/M — 
(1 - (£:7L)2/3)-3/2^ corresponding to /?oo = {£'^/Lf'^. For fixed M and £ the minimum 
energy loss is reached at /S^m- Assumption (3oo — Pmin, then gives 

{l + iPooTooaooff"^^^ (All) 
which in the strongly relativistic limit gives Michel solution 

Too = (A12) 

We can also relate the terminal magnetization a^o to the magnetization at the source - 
more specifically to magnetization at the sonic point o"/: 

ScToo if (Too -C 1 
= IcToo if (Too > 1 



Thus, we always have (Too < ct/j but they remain of the same order of magnitude: the mag- 
netization of the fiow changes only slightly as the fiow propagates away from the launching 
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point to infinity. The reason for constant a in the supersonic regime is that both in the case 
p ^ p (hnear acceleration stage) and p <^ p (coasting stage) the plasma and the magnetic 
field energy densities in the flow change with the same radial dependence (~ and ~ 
correspondingly) . 

For arbitrary flow parameters the evolution equations are integrated numerically (Fig. 
3). Given the evolution of the flow and the relation for local a we can find the evolution of 
the magnetization parameter (Fig. 3.b). 
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surfaces 

Fig. 1. — Cartoon of an expanding, coUimated flux rope. Twisting of footpoints of a flux 
tube leads to electric current flow along the loop and results in its expansion. Beyond some 
twist angle dynamic instability leads to formation of dissipative current sheets at the leading 
edge of the loop, between the footpoints and, possibly in the bulk. Reconnection at the 
leading edge allows the flux tube to break out of the magnetosphere. CoUimating effect of 
the tied footpoints later lead to broadly coUimated outflow. 
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Fig. 2. — Doppler factor V — l/(r(l — (5cos9) (solid line) and apparent transverse velocity 
(dashed line) for a relativistic shock expanding with Lorentz factor F = 3 and moving at 
6ob — 135° with bulk Lorentz factor V^uiy. = 2 as a function of the angle 9 between the 
explosion direction and the emission point in the rest frame of the explosion for emission 
points located in the plane containing direction of bulk motion and line of sight. 
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Fig. 3. — Evolution of strongly magnetized and hot flows, (a) Four-velocities of flows are 
given as functions of r/r/ for L/M — 100 and different values of the parameter ^. Flows 
start at r = r/ with (5 — (5f] supersonic flows first accelerate as ~ r, reaching a terminal 
value given by the larger root of eq. (AlO), while subsonic decelerate initially as /3 ~ r"^ 
reaching asymptotic value given by the smaller root of eq. (AlO). (b) Magnetization param- 
eter (7. For supersonic flows (lower branch) the magnetization remains constant, reaching 

Coo — ~ ^z) asr— >oo. Subsonic flows become strongly magnetized as they expand 
(upper branch); the magnetization parameter increases a ~ r^^^. 



